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Recent advances in submicron technology have made it possible to fabricate ballistic 
quantum dots in which the electron mean free path exceeds the system size. For dots with 
irregular shapes, the dynamics is largely chaotic |1| and its signatures are observed in the 
transport properties of the device 0. In dots that are weakly coupled to external leads 
via tunnel barriers, the resonances are isolated, and at low temperatures the conductance 
is dominated by the resonance that is closest to the Fermi energy in the leads. Such closed 
microstructures thus offer a unique opportunity to probe the chaotic signatures of individual 
wavef unctions. Since the charging energy of the dot is large compared with the mean-level 
spacing, the conductance exhibits a series of almost equally spaced peaks versus the gate 
voltage (Coulomb blockade oscillations). The peaks show order-of- magnitude fluctuations 
which are explained by a statistical theory 0,|J. The distributions of the conductance peak 
heights were measured for both conserved and broken time- reversal symmetry |5]||, and were 
found to agree with theory. Using the supersymmetry technique for disordered systems f?||| 
and random matrix theory (RMT) ||, closed expressions for width and conductance peaks 
distributions were derived for any number of possibly correlated channels for both orthogonal 
and unitary symmetries. 

In this letter we derive exact expressions for the universal width and conductance peak 
distributions in the crossover regime between conserved and broken time-reversal symmetry 
for any number of possibly correlated and/or inequivalent channels. An approximate expres- 
sion for the distribution of the wavefunction intensity at a fixed spatial point was obtained 



in Ref. ||10|| , and an exact expression was derived in Ref. [|TTJ using supersymmetry. The 
averaged intensity over the complete spectrum of the random matrix was derived in Ref. 
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At low temperatures T -C kT -C A, the conductance peak amplitude G is given by (13 



G = 4 where 
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Here T 1 ^ is the width of a resonance to decay into the left (right) lead and Y is the total 
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average width, g is dimensionless and temperature- independent. In general we assume that 
the left (right) lead has A^ r ) open channels such that T 1 ^ = I] c |7^ r - ) | 2 , where 7^ r ) is the 
partial amplitude to decay into channel c on the left (right). Using R- matrix theory [[| or 
resonance theory ||14|| , the partial width amplitude can be expressed in terms of the resonance 
wavefunction \l/ and the channel wavefunction $ c at the lead-dot interface. When each lead 
is modelled in terms of several point-like contacts r c then 7 C oc ^(r c ) . By expanding a 
resonance wavefunction in a complete basis of solutions p^if) inside the dot at a fixed energy 
(\&(r) = Vv/V( r ))> the partial width j c can be expressed as a scalar product fL4| of the 



vectors that represent the channel and the resonance function 7 C = (4> c \tp) = X^^^/V 

In the crossover regime of breaking time-reversal symmetry, the resonance wavefunc- 
tions are assumed to have statistical properties which are described by the corresponding 
eigenvectors if) of a random matrix ensemble that interpolates between the GOE and GUE 
I 

H = S + iaA . (2) 

S and A are, respectively, symmetric and antisymmetric real matrices of dimension N which 
are uncorrelated and chosen from gaussian ensembles of variance a 2 . The proper transition 
parameter that describes the crossover from GOE to GUE is given by the ratio of the rms 



of a typical symmetry-breaking matrix element to the mean-level spacing A |L5 



aa ayN , , 

where we have used the mean level density at the middle of the spectrum. The complete 
breaking of time- reversal symmetry occurs for A ~ 1. Since the transition parameter depends 
on the level density, we shall study the ensemble's statistics only around the middle of 
the spectrum. The spectral properties of the transition ensemble (0) were derived in a 



closed form [111. Similar spectral correlators were also derived for a single electron in a 



disordered medium using the supersymmetry method |T7|]. However, less is known about the 
eigenvector statistics and, unlike in the GOE and GUE limits, the eigenvalue and eigenvector 
distributions do not factorize. 



An eigenvector component is a complex number, ip^ = ip^R + iip^i-, and can be viewed 
as a two-dimensional vector in the complex plane. Since the eigenvector if} is determined 
only up to a phase e l9 , ip^ is determined only up to a rotation by an angle 9. This angle is 
uniquely determined by rotating to a principal frame in which |p"5| . |T0f 

N 

E Vv*Vv = o ; E E = ^ 2 ( 4 ) 

with < r < 1. The ratio r 2 is invariant under a rotation and is thus independent of 
the frame. The right inset of Fig. 1 shows the components of a typical eigenvector in the 
complex plane. Its general shape is that of an ellipsoid whose semi-axes define the principal 
axes. In the following, ip^R and ip^i will denote exclusively the components of ip in this 
principal frame. 

Let us consider eigenvectors with a fixed value of r. Under an orthogonal transformation 
(in the iV-dimensional space) the real and imaginary parts of ip do not mix so that a 
principal frame remains principal and r 2 is invariant. Since the probability distribution 
of the ensemble (0) is also invariant under an orthogonal transformation, the conditional 
probability of finding an eigenvector t/j given its r-value is 

The conditional distribution P(X\r) of any quantity X which is a function of the eigenvector 

■0 can be calculated from (||). The actual distribution P\(X) at a given value A of the 

transition parameter can then be calculated in terms of the distribution P\{r) of the quantity 

r (whose explicit form is given below in Eq. (||)) 

i 

P X (X) = J P x (r)P(X\r) = (P(X\r)) . (6) 

Q 

The conditional joint distribution of the partial width amplitudes 7 = (71,72, • • • ,7a) 
for A real channels 4> c can be calculated exactly from (|5|) and 7 C = (<fi c \ift). In the limit 
iV —>■ 00 we obtain a Gaussian distribution in 7 C # and 7 C / 

P( 7 |r) ex exp (l^-^M- 1 ^ + l^^M' 1 ^] , (7) 
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where 7 C = j c r + ij c i = {4> c \ipn) + i((ft c \ipi) is a principal frame decomposition, and 7 T 
denotes the transpose of the column vector 7. The matrix M is the channel correlation 
matrix M cc / = 7*7c' = (<t> c \4> c i) ■ To calculate M for a chaotic system, we note that the 



relation V'uVy = iV -1 ^/ is valid for the transition ensemble irrespective of the value of 
A. We can therefore repeat the derivation of Ref. || to find for M an expression identical 
to the one obtained in either the orthogonal or unitary limits. For point-like contacts 
M cc i oc Jo(k\r c — r^l), independently of the crossover magnetic field. 

As a special case of (0) we obtain the conditional distribution of a single component ip^ 
of an eigenvector by choosing a single channel along the /x-th axis. At a fixed r, ip^R and 



ip^i are independent Gaussian variables with ip^j/ip^R = r 2 . The full distributions of the 
real and imaginary parts of the wavefunction amplitude are obtained by a weighted average 
of the conditional distributions according to Eq. (|6]) and are not statistically independent. 
Fig. |l| shows these distributions for A = 0.1. A Gaussian distribution (dashed line) is a good 
approximation for the real part of the wavefunction amplitude, but not for the imaginary 
part. 

The distribution of the width r = |7| 2 = 7^ + 7I for a one-channel lead is found to be 

P x (t) = (a + e- a + f / (a+a_f)) , (8) 

where T = T/T, a± = (r -1 ± r)/2 and 1$ is the modified Bessel function of order zero. An 
exact expression for the distribution of the wavefunction intensity (at a fixed spatial point) 



was derived in Ref. [11] through the supersymmetry method. By comparing our result 



with the result in [11], we find that the r-distribution is given by 



P A (r) = vr 2 A 2 (l/r 3 - r)e-^ A2(r " 1/r)2 {^(A) + [(r + l/r) 2 /4 - 1/2tt 2 A 2 ][1 - ^(A)]} , (9) 



where 0i(A) = j e 27I " 2a2 ( 1 f2 ^dt. RMT simulations confirm that Eq. (P) is indeed the r- 


distribution (see left inset in Fig. [1]). Eqs. (0), (|^) and (^) provide a closed analytic 
expression for the joint partial amplitudes distribution for a lead with any number of possibly 
correlated and inequivalent channels. 



To derive a closed expression for the total width distribution for multi-channels leads, we 
note that the conditional distribution P("y|r) in Eq. (|7|) is identical to a GOE distribution 
for 2A channels with partial amplitudes 7 C #, and a correlation matrix M. composed of 
four A x A blocks 



M 



' -^—M » 



l+r 2 



V TT^Mj 



(10) 
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We can therefore use the known distributions from the GOE limit [|9|. The 2A eigenvalues 
of M. are given by {uj 2 } = {j^w 2 , j^w 2 } where w 2 are the A eigenvalues of M. Sorting 
the inverse eigenvalues of Ai in ascending order uj^ 2 < u^ 2 < . . ., we have 
1 /■' , e~ Tr 

i2m-l( _ = T 

Similarly we can use the closed expression for the GOE conductance distribution || 
to obtain P(g\r) analytically and then use (|6|) to find P\(g) = (P(g\r)). In this calcu- 
lation of P(g\r) (where Eq. (|I|) is exploited), we take advantage of the statistical inde- 
pendence of the conditional distributions of the total widths in the left and right leads, i.e. 
P(r , r r |r) = P(T l \r)P(r r \r). We note however that the widths themselves are not indepen- 
dent since (P(T l \r)P(T r \r)) ^ {P(T l \r))(P(T r \r)). The distant correlation of wavefunctions 



in the crossover regime |1| follows from this relation. The universal conductance peaks dis- 
tributions P\{g) for one-channel symmetric leads are shown in Fig. ||] for several values of A 
in the crossover regime. A good approximation (see left inset in Fig. |2|) is P\{g) ~ P(g\ro), 
where r = ?"o(A) (see right inset in Fig. 0) is determined by finding the best fit to the exact 
P\(r) for a single channel. r$ can also be estimated by (r) (solid line in the right inset in 
Fig- 0). 

To test our RMT predictions for the crossover conductance distributions, we used the 
conformal billiard ]0| (threaded by an Aharonov-Bohm flux line $) whose shape is deter- 
mined by the image of the unit circle in the complex z-plane under the conformal mapping 



w(z) — (z + bz 2 + be lS z 3 ) / y/T+Eb 2 . We collected statistics from several uncorrelated and 
fully chaotic billiards with 5 = it/ 2 and various values of b. Semiclassical considerations 
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for weak fields p0[ lead to a linear relation A = $/$ cr . $ cr is the crossover flux given by 
^cr/^o = ( 7r /4)(agA/')~ 1//4 , where A/" is the number of electrons within the dot and a g is a 
geometrical factor. Except for a constant factor, this crossover field has a similar expression 



to that of the correlation field ||14|| . Rather than calculating the geometrical factor semiclas- 
sically, we determined the exact relation between A and $/$o by fitting the Dyson-Mehta 
A3 spectral statistics to its known analytic form [jnyHJ. The results for the lowest 300 eigen- 
values are shown in the top inset in Fig. For each lead we chose a sequence of A equally 
spaced points on the billiard boundary with spacing |Ar|. To test our expression for the cor- 
relation matrix M, we have calculated the averaged wavefunction amplitude correlations and 
found them to be in good agreement with J (k\Ar\) independent of the flux (see the inset 
to the middle panel of Fig. 0). Fig. |] compares the conductance distributions in the con- 
formal billiard with $/$o — 0.04 (histograms) to the theoretical predictions (using A = 0.16 
from the inset) for three cases: one-channel leads, two-channel leads with fc|Ar| = 0.5 and 
four-channel leads with fc|Ar| = 2.4. In each case we also show the limiting distributions for 
conserved and fully broken time-reversal symmetry. The billiard calculations are in good 
agreement with the intermediate distributions. 

In conclusion, we have derived in closed form the width and conductance peak distribu- 
tions in a chaotic quantum dot for leads with any number of possibly correlated channels in 
the crossover regime from orthogonal to unitary symmetry. The distributions depend only 
on the symmetry-breaking transition parameter and the channel correlation matrix M in 
each lead. This work was supported in part by the DOE, Grant No. DE-FG02-91ER40608 
and by the NSF, Grant No. PHY94-07194. NDW was supported by the European Hu- 
man Capital and Mobility Programme and by NSERC Canada. YA acknowledges useful 
discussions with S. Tomsovic. 
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FIGURES 

FIG. 1. The distributions P{^hr) (wider) and P(ipui) (narrower) of the real and imaginary 
parts of an eigenvector component in the principal frame for A = 0.1. Histograms: simulations of 
the RMT ensemble (||); dashed lines: a Gaussian approximation; solid lines: the exact analytic 
results. The right inset shows the components of a typical vector in the complex amplitude plane. 
The left inset is the r-distribution P\(r) for A = 0.1, where the solid line is the analytic result (|9|) 

and the histogram is from RMT simulations. 

FIG. 2. Conductance peak distributions P\(g) vs. log g in the crossover from conserved to bro- 
ken time-reversal symmetry for one-channel symmetric leads: A = (GOE, dashed), 0.1,0.25,0.5 
(solid lines) and A> 1 (GUE, dot-dashed). Curve maxima increase with A. The left inset compares 
the Gaussian approximation for P{g) (dashed) with the exact result (solid) for A = 0.1. Shown in 

the right inset is ro(A) (circles) and (r)\ versus A (solid line). 

FIG. 3. Conductance peak distributions P(g) in the conformal billiard (histogram) for A 

point-contact symmetric leads and flux of ^/^o = 0.04 for several values of A and fc|Ar|. In 

each case we show the analytic predictions (solid lines) as well as the GOE (dashed) and GUE 

(dot-dashed) limits. The top inset describes the transition parameter A as a function of magnetic 

flux ^/^o- The inset in the middle panel shows the spatial correlations of the eigenfunctions for 

= 0.02, 0.06 and 0.10 (diamonds, pluses and x's, respectively). The solid line is the theoretical 

prediction. 
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